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perturbative results for the splitting functions, we present our scheme with the complete contribution 
from the double logarithms, being the largest soft gluon logarithms. We show that the resulting 
approximation is more complete than the Modified beading bogarithm Approximation even with 
the fixed order contribution calculated to leading order only, and find, after using it to fit quark and 
gluon fragmentation functions to experimental data, that this approximation in our scheme gives 
a good description of the data from the largest x p values to the peak region in £ = ln(l/x p ), in 
contrast to other approximations. In addition, we develop a treatment of hadron mass effects which 
gives additional improvements at large £. 
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I. INTRODUCTION 

The current description of single hadron inclusive pro- 
duction processes within the parton model of perturba- 
tive QCD (pQCD) is provided by fragmentation func- 
tions (FFs) D^(x, Q 2 ), each of which corresponds at low- 
est order to the probability for the parton a produced at 
short distance 1/Q to form a jet that includes the hadron 
h carrying a fraction x of the longitudinal momentum of 
a. Different theoretical schemes have been derived de- 
pending on the kinematic region of x: Fixed order (FO) 
calculations at intermediate and large x and resumma- 
tion to all orders of soft gluon logarithms (SGLs) at small 
x. What is needed is a single formalism valid over the 
union of all ranges that the different pQCD approaches 
allow. This unification must be consistent, i.e. it must 
agree with each approach in the set, when the expansion 
of that approach is used, up to the order being consid- 
ered. 

Much progress has been made 0, m determining 
fragmentation functions (FFs) at large and intermedi- 
ate momentum fraction x using FO Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) evolution [3( to next- 
to-leading order (NLO) [j]. However, determination of 
FFs at small x is performed independently, since the cal- 
culation of the evolution requires a different approach 
because of soft gluon logarithms (SGLs). The complete 
resummed contribution from the largest SGLs, being 
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the double logarithms (DLs), of the splitting functions 
is obtained from the Double Logarithmic Approxima- 
tion (DLA) d 13, while some information on the con- 
tribution from the next largest class of SGLs, the single 
logarithms (SLs), in the splitting functions is obtained 
from the Modified Leading Logarithm Approximation 
(MLLA) [10,I3- (The complete evolution is obtained by 
using the approximation that the quark FFs are identi- 
cal, and equal to the gluon FF multiplied by a calculable 
constant.) Provided all incomplete higher order terms 
are not allowed to become too large, the MLLA evolu- 
tion can describe small x data very well However, 
although the cross section over a large range of x can be 
described by pQCD in general, two different approaches, 
for which no matching conditions exist, are required to 
relate one to the other. In addition, the range of x over 
which both approaches are valid is not clear. 



In this paper, we are concerned with the unification 
of DGLAP evolution and the resummation of SGLs in 
pQCD. In Sec. [HI we derive a simple and consistent 
scheme which reproduces both approaches to any desired 
order. In Sec. IIIII we use the DLA to obtain the complete 
DL contribution to the evolution within this scheme. 
These two sections justify and extend the scheme intro- 
duced in Ref. [l(|. In Sec. II VI we compare this approx- 
imation with experimental data. Finally, in Sec. we 
present our conclusions. For the reader's convenience, we 
list the explicit results for the leading order (LO) split- 
ting functions in the FO approach in Appendix^ as well 
as references for the explicit results for the NLO ones. 
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II. SGL RESUMMATION IN DGLAP 
EVOLUTION 

In this section, we discuss DGLAP evolution and its 
SGLs in the FO approach, and then we give a formal 
definition of our general scheme in which the SGLs in 
this approach are resummed. 

The DGLAP equation reads 

^—D(x,Q 2 ) = j ld Vp( z ,a s m)D( X -,QA, (1) 
dlnQ z J x y \V J 

where, for brevity, we omit hadron and parton labels. D 
is a vector containing the gluon FF D g and the quark and 
antiquark FFs D q and Dq respectively, in linear combi- 
nations according to the choice of basis, and P is the ma- 
trix of the splitting functions. We define a s = a s /(27r), 
whose Q 2 dependence is determined by the QCD (3 func- 
tion /3(a s (Q 2 )), through the Callan-Symanzik equation 

^A_a s {Q 2 )^P{a s {Q 2 )). (2) 

The (3 function can be calculated in perturbation theory, 
where it takes the form 

oo 

/?(a s ) = - X>n-2<- (3) 

n=2 

Choosing a factorization scheme in which P is explicitly 
independent of quark masses and then exploiting the re- 
sulting SU(n/) symmetry for n/ quark flavours, and also 
exploiting the charge conjugation invariance in P, leads 
to the simplest basis for writing Eq. Q, which consists 
of the combinations (i) D = D~ , where D~ = D q — Dq 
is the valence quark FF, (ii) D = Dns> where Pats is a 
non-singlet quark FF, i.e. any linear combination of the 
FFs .D+ = D q + Dq which vanishes when they are all 
equal, and (iii) D = (Ds,D g ), where 

n f 

D ^ = -J2 D t ( 4 ) 

is the singlet quark FF. The n./ — 1 non-singlets must 
be chosen such that together with the singlet a linearly 
independent basis for the quarks is formed. 

We will often work in Mellin space, where a function 
f(x) becomes 

/H = / dxx"f(x), (5) 
Jo 

where u) is any integer greater than those values of lu 
for which f(u>) is non-singular, since the convolution in 
x space in Eq. JJ), and in equations later in this paper, 
becomes the simple product 

-r^Diu, Q 2 ) = P(u, a s (Q 2 ))D(uj, Q 2 ). (6) 



After performing the desired analytic operations, Mellin 
space results can be transformed back to x space by an- 
alytically continuing f(ui) to complex ui and then using 
the inversion formula 

/(*) = ^ / dux-"- 1 /^), (7) 

where C is a contour in Mellin space from Im(cj) = — oo 
to Im(ui) — oo, which passes to the right of all poles in 

/H- 

Without knowledge of P, Eq. (or Eq. ©) provides 
no constraint on D(x, Q 2 ) over the ranges < x < 1 and 
< Q 2 < oo. For a given u>, specification of P(u>, Q 2 ) for 
all Q 2 fixes the the Q 2 evolution of D(oj, Q 2 )- This means 
that, given D(uj, Qq) at some specified Qq, D(uj, Q 2 ) can 
be calculated for all Q 2 . This evolution is usually calcu- 
lated explicitly for the evolution matrix E, defined for all 
Q 2 and Qq by 

D(w, Q 2 ) = E (w, a s (Q 2 ), a s (Q 2 )) D(u, Q 2 ), (8) 

and for all a s , ai and clq by 

E(uj,a s ,a s ) = I, ^ 
E(w,a B ,ai)E(u,ai,a ) = E(uj,a s ,a ), 

(I is the unit matrix) with no loss of generality for the 
functional form of D(oj, Q 2 ). Then, with some additional 
definitions for E that result in no loss of generality, E 
becomes fully constrained in terms of P by invoking Eq. 
©, 

P{uj,a s {Q )) = d\ n Q2 E \u,a s {Q ),o ). 

(10) 

The boundary conditions in Eq. JUJ) can be used to verify 
that the right hand side of Eq. 110|) is independent of ao- 
The factorization theorem 1 1 _1| states that P is an in- 
variant with respect to the hadron being observed, and 
furthermore that the series for P{x,a s ) in a s keeping x 
fixed, 

oo 

P(x,a s ) = Y,<P (n ~ 1} (x)i (11) 
n=i 

can be calculated from perturbation theory even when 
any quark masses go to zero or infinity. Equation 
truncated at some chosen (finite) n is known as the FO 
approach, and is not valid at small x due to the pres- 
ence of terms which in the limit x — > behave like 
(a™ / x) m 2n_m_1 x f or m = 1, In— 1. Such logarithms 
are called SGLs, and m labels their class. As x decreases, 
these unresummed SGLs will spoil the convergence of the 

— 1/2 

FO series for P(x,a s ) once ln(l/x) = 0(a s ). Conse- 
quently the evolution of D(x, Q 2 ) will not be valid here, 
since the whole range x < y < 1 contributes in Eq. 
Therefore, the FO approach is only a good approximation 
for sufficiently large x. 
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SGLs are defined to be all those terms of the form 
m only, where m = 1, 2n and labels the class 



of the SGL, in the expansion about u> 
transform of Eq. , 



of the Mellin 



P{u,a s ) 



(12) 



We will consider SGLs of the type to — 2n later. For 
to = 1, 2n — 1, this definition agrees with the form of 
the SGLs in x space given above, since 



1 

UJP 



(-l)f 



dxx u 



(13) 



for Re(w) > and p > 1. Such terms spoil the conver- 
gence of the series in Eq. Ijl2|l as w — » 0. What we require 
is an alternative scheme for the evolution which will be 
valid for large and small to. The inverse Mellin transform 
of this evolution should then be valid for large and small 
x. For this purpose, we propose the following general 
scheme, which we call the SGL+FO scheme. Firstly, P 
is written in the form 



p _ pFO + pSGL 



(14) 



where P SGL contains only and all the SGLs in P, so that 
P FO is completely free of SGLs. Secondly, by summing 
all SGLs in each class to, P SGL (w, a s ) is resummed in the 
form 



(15) 



and truncated for some finite m. The functions g m {x) 
in Eq. (|15f) are not Taylor series in either x or any func- 
tion thereof. Note that, apart from the condition that 
the series must start at m — 1, which follows from the 
definition of SGLs above, Eq. I|15|) is just the general re- 
sult of expanding a function of a s and U) in a s /u keeping 
a s /uj 2 fixed. Thirdly, the remaining FO contribution to 
P, P (lu, a s ), is expanded in a s keeping oj fixed, 



P FO ( W ,a s ) = X:< J P FO(n - 1) M ! 



(16) 



n=l 



and truncated for some finite n. P FO (ui, a s ) can be ob- 
tained by subtracting all SGLs from the series for P on 
the right hand side (RHS) of Eq. 1|12|) . Since all classes 
m < 2n are included in p SGL j i.e. since P SGL contains all 
terms of the form a™ /ui p for p = 0, 2n — 1, P fo (uj, a s ) 
is zero when to = 0, since 

pFO(n) ( Q j = o Fourthly and fi- 
nally, the result for P(lo, a s ) is inverse Mellin transformed 
to obtain P(x,a s ), and then Eq. is solved exactly 
(which can be done numerically). 

From (incomplete) calculations of Eq. I|15|l up to the 
class m — 2 &}, P is believed to be finite at uj = 0, and 
in particular to be a series in Ja^ with finite coefficients, 



beginning at O(^fal). This means that each term of the 
form a™/tu p for the types p > 1 in the expansion of Eq. 
(I12JI about u> = should be included in the resummed 
term of class to = 2n — p. However, terms of the type p = 
(m = 2n), which are included in our definition of SGLs, 
are non-singular and may therefore be left unresummed. 
Thus we define a second general scheme, which is the 
same as the SGL+FO scheme defined above, except that 
we separate 

pSGL in Eq ^ intQ 



3SGL 



_ pSGL 



P, 



SGL 



(17) 



and expand only P^> ^ as in Eq. I|15[l , while Pp3o > which 
is independent of w, is expanded as a series in a s . We 
shall call this the SGL+FO+F05 scheme, where "+FO<5" 
means that the p = terms, which are each proportional 
to 5(1 — x) in x space, are left as a FO series in a s . 

To summarize our SGL+FO (+FO<5) scheme, we resum 
SGLs for which m = 1, 2n (m = 1, 2n — 1) in the 
form of Eq. i|15|) , and treat all remaining terms as in the 
FO approach. 

In the phase space region for a s < 1 and x above val- 

— 1/2 

ues for which ln(l/x) = 0(a s ), the following x space 
results indicate that the SGL+FO (+FO<5) scheme gives 
a good description of the evolution. Firstly, Pp>i, ob- 
tained from the inverse Mellin transform of Eq. (|15[) with- 
out terms of the type p — 0, can be written as 

oo 

P*§ l(x, a.) = ^ lna; ) m f™ ( a * ln ' x ) ■ ( 18 ) 

rn — l 

Equation (|18|l can also be obtained by summing the SGLs 
in x space for each to. Since a s hix is always small, the 
series in Eq. l|18fl is a valid approximation when x is small. 
On the other hand, as x — *■ 1 the SGLs for the types p > 1 
all vanish, and therefore so does each term in the series 
in Eq. I|18|) . Secondly, the full contribution from the type 
p = terms is just 



P^(x,a s ) =5{l-x)Y J C n t 



(19) 



n=l 



Thirdly and finally, the expansion of P FO (x,a s ) in a s , 
i.e. the inverse Mellin transform of Eq. I|16|). 

oo 

P FO (x,a s ) = J2<P F ° in - 1) (x), (20) 

71=1 

converges for all x. 

III. DLA IMPROVED DGLAP EVOLUTION 

In this section, we summarize the DLA and its re- 
lation to the DGLAP equation, and use this under- 
standing to derive a modified form of DGLAP evolu- 
tion in which all DLs, the to = 1 class of SGLs, are 
included and resummed in a manner consistent with the 
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SGL+F0(+F0<5) scheme denned in Sec-EJ The DL con- 
tribution to the evolution obeys the DLA equation, which 
can be obtained from the DLA master equation for the 
quark and gluon generating functionals given in Ref. 
and is given by 

Jx v v \y 



x V V 



a s (Q 2 )D[- - 2 
V 



(21) 



(The values of various factors are given in Appendix |A"1) 

2 d 

The result y dln « 2 f{Q 2 ) = f{y 2 Q 2 ) has been used to 
obtain the second line in Eq. (|21|l . Explicitly A = for 
the DL evolving parts of the components D = D~ and 
D = Dns, while 



.4 



r\ 2CV 
U C A 

1 



(22) 



for the component D = (D^,D g ). Note that A is a 
projection operator, i.e. it obeys 



A 2 = A. 



(23) 



In Mellin space, Eq. I|21|l . with certain boundary condi- 
tions to be discussed later, is equivalent to Eq. JQl with 
P replaced by the m = 1 term in Eq. (|15|l , 



-.91 



(24) 



up to incomplete higher order terms in the remaining 
FO contribution and up to incomplete SGLs of classes 
for which m > 2. 

The remaining part of the evolution (i.e. all the FO 
terms and remaining SGLs) can be included in Eq. I|21|) 
by writing it in the form 



d 



d\nQ- 



;D(x,Q 2 ) 



a s (Q 2 )D ( -,Q 2 
V 



y y 

" 1 ^p( y ,a s (g 2 ))^f- 
y \y 



(25) 



P(x,a s ), which must be free of DLs, is constrained in 
terms of P since Eq. Ij25(l must be equivalent to Eq. Q 
(for D = D~ and D = D^g, one obtains the trivial 
result P = P). In general, this can be done explicitly by 
expanding the operator in Eq. lll'ol) in the form 



2 



exp 



21ny 



dhiQ 2 



y> (21ny)" ( d 

71=0 



dlnQ 2 

(26) 



and then repeatedly applying the evolution equations, 
Eqs. Q and ©, to the (ot^)" [a s (Q 2 )D (f , Q 2 
operations in Eq. (|25|l . For example, to 0(a 2 ) one finds 



P(x,a s ) = P(x,a s )-2C A A 
1,1 dy 2a s j;m | 



a s Anx 

— + 2(3(a s ) 

x x 



P(.V,a a ) 



(27) 



+ 0(a 3 s ). 



In the square brackets on the RHS of Eq. (|27(l . only 
the first term contributes to the 0(a s ) (LO) part of P, 
while the second and third term contribute to the 0(a 2 ) 
part. To this accuracy, the third term is calculated with 
P(x,a s ) = a s pW(x), which can be found in the litera- 
ture (see Appendix . From Eq. H27|) we observe that 
P up to 0(a 2 ) (NLO) is free of DLs, by taking the DLs 
in P(x,a s ) for the first term on the RHS of Eq. f2T|l up 
to 0(a 2 ), 



P UL (x ias ) 



2CA—a s 

x 



4c2 Aln 2 ^ 2 



ai + 0(ai). (28) 



In practice, it is numerically very difficult to solve Eq. 
B25JI (and Eq. I|21(l ). in particular when x < Qq/Q, which 
requires calculating the FFs for Q < Qo- This prob- 
lem corresponds simply to the fact that in addition to 
knowing D(x,Qq), which is required to solve Eq. JQl, 

,, d „.i D(x, Q 2 )\ _ „ must also be known in order to 

^ .Z— 'Q=Qo __ 
solve Eq. I|25() . Even more seriously, Eq. I|25[) cannot 

be solved numerically at all when x < Aqcd/Q, due to 
the Landau pole in a s (y 2 Q 2 ). 

Instead, we will examine what constraint Eq. i|25[) pro- 
vides for P, since then the evolution can be performed 
using Eq. Q), which is numerically easily solved, with- 
out requiring explicit use of D(x,Q 2 ) at scales less than 
Qo at any x. For this purpose, we work in Mellin space, 
where Eq. I)25|l becomes 



uj + 2- 



dlnQ 2 J rflnQ 2 



uj + 2 



D{lo,Q 2 ) = 2C a cl s {Q 2 )AD(lo,Q 2 ) 



P(Lu,a s {Q 2 ))D(io,Q 2 ). 



dlnQ 2 



(29) 



After substituting Eq. © into Eq. 129|) and dividing out 
the overall factor of D(ui,Q 2 ), we obtain the following 
constraint on P: 



LU + 2 



d 



d\nQ- 



(P - P) + 2 (P - P) P - 2C A a s A = 0. 

(30) 

Equations (|25|l . I|29(l and H3()|l are exactly equiva- 
lent, however Eq. I|3U|I shows most clearly that, for all 
Q 2 , specifying P(u,a s (Q 2 )) will completely constrain 
P(uj,a s (Q 2 )) once P(to, a s (Qo)) is chosen. However, the 
only information we have for P is that it is free of DLs, 
which means that an explicit constraint can be obtained 
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only for P DL , which we now do. We first make the re- 
placement 

P = P + P DL (31) 
in Eq. (|30|l . where, comparing with Eq. l[Tl)l. 

p = pFO + pSGL _ pDL (32) 

Then we expand Eq. H3()|l as a series in a s /uj keeping 
a s /oj 2 fixed and extract the first, 0((a s /u>) 2 ), term to 
find that the constraint on P DL is exactly 



2(P DL ) 2 + loP dl - 2C A a s A = 0. 



(33) 



Equation (|33|l gives two solutions for each component of 
P. Since P is never larger than a 2x2 matrix in the 
basis consisting of singlet, gluon, non-singlet and valence 
quark FFs, there are four solutions which read 



P^(Lj,a s ) = S±(u,a s )A, 
P^(oj,a s ) = -^I-S ± (u;,a s )A, 



(34) 



where 



5±(w, a s ) = \ (-w ± V^ 2 + 16C A a s ) . (35) 

For the evolution matrix E defined in Eq. (JSJ), these so- 
lutions for P DL correspond respectively to 



(36) 



(38) 



(39) 



where we have used Eq. (|23|) . and defined 

R±{cj,a si a a ) = / da — , 

/3(a) 



P(w, a S) ao) 



a; f aa da 



2A /9(a) 



Here, Eq. (J2J) has been used to transform the In Q 2 inte- 
grals into integrals over a s . The general solution to the 
DL part of Eq. Q is, finally, Eq. © with 

P(w,a s ,a ) = yj^t,j(a;, a s , a )ki,j(w, ap), (40) 



where i = 1 , 2 and j — ±. To ensure that E is normalized 
as in the first line in Eq. (jSJ), the matrices fc^-t obey 



^ • fcjj (a>, ao) = I. From Eq. (fTTTf) . the most general 
splitting function is then 



P DL (oJ,a s ) = 



^T^' o*)- 5 ^ (^' a s , a )fci,j (w, ao) 



E it j(u}, a s , a )ki tj (u), a ) 



(41) 



The only solution for P DL (w, a s ) which is consistent with 
the DLs in the known result for P in the FO approach 
of Eq. l(TT|l to 0(a 2 ) (see Appendix is that for which 

= / and all the other kij are zero, i.e. P DL = 
or more explicitly, 



A 



P DL (w, a s ) = — (-w + \Juj 2 + 16C A a s 



(42) 



since, in the component D = (D^,D g ), the expansion of 
the result in Eq. l|4^|) in a s to 0(a 2 ) keeping u fixed gives 



P DL (^,a s ) 



as icv _ a 2 







2Ca 



,28Cl 



O(aJ), (43) 



i.e. the Mellin transform of Eq. I)28fl . while in the com- 
ponents D = D~ and D = D NS , P DL = 0. The other 
possibilities implied by Eq. Q41JI do not give these re- 
sults and/or cannot be expanded in a s , i.e. they contain 
non-perturbative terms. Equation l|42|l agrees with the 
results of Ref. 0, which are derived using the conven- 
tional renormalization group approach, and with the re- 
sults from the generating functional technique of Ref. • 
Thus, the explicit evolution is given by Eq. JSBJl with the 
upper sign. Writing R + = R and returning to Eq. JSJ, 
the evolution obeys 



D(u, Q 2 ) = I + [e R ^-^ Q *).».»8)) - lj A D(w, Q%). 

(44) 

For all DLs in the evolution, Eq. (glj) solves Eq. (|2*T)l . 

A completely explicit form for Eq. I|44(l can be obtained 
in the case that Eq. (JSJ is taken to 0(a 2 ) only. Equation 
© then implies that a s (Q 2 ) = l/(P \n(Q 2 /A 2 QCI) )), for 
which Eq. (|38|l (with the upper sign) reads 



R(uj : a s ,a ) -- 
2C A 



4ft 



- (-U3 + \fb. 

oa s V 



16CAO; 



In 



(45) 



(a s <-> ao). 



w/3 "* -w + \/lo 2 + 16CUa s 

Since Pgt = P^ = according to Eq. g2J|, it follows 
from Eq. © when P = P DL that, for small u>, 



d\nQ 



;D^(u,Q 2 ) 



Pg g L (a,,a s (Q 2 )) d 
P^(u;,a s (Q 2 ))d\nQ 2 



D g (co,Q 2 ). 
(46) 
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With the results for P|? g L and P^ L in Eq. (JUJ, integrating 
Eq. l(4l)fl over In Q 2 gives 



2C? 



(47) 



The constant of integration has been neglected in Eq. 
(P7Jl . which is valid for large Q. Equation I^Tfl (again up 
to an additional constant) can also be derived from Eq. 

Since the non-singlet and valence quark splitting func- 
tions are free of DLs, the derivatives of the non-singlct 
and valence quark FFs with respect to In Q 2 may be ne- 
glected at small ui. Again, integrating such results over 
In Q 2 and neglecting the constants of integration implies 
that the non-singlet and valence quark FFs vanish. In 
this case Eq. I|47JI becomes 



D„ 



(48) 



reducing the number of FFs required for the cross section 
to one, D g . Such a low x approximation is often used in 
DLA or MLLA analyses of data. However, since we want 
a complete formalism suitable for both small and large x, 
we will only use Eq. (|48|l to partially constrain our choice 
of parameterization at low x in the next section. 

Complete information on the SL contribution to P SGL , 
given in the notation in Eq. I|15l) by 



(49) 



cannot be obtained from Eq. (|30|l . since the full SL con- 
tribution to P is not known. However, its SL at 0(a s ), 
which according to Eq. l|2*7|l is equal to the SL in P at 
0(a s ), given according to Appendix lAl bv 



o 



-3C F 



l T R n f 



(50) 



is a type p = term (see Sec.|Hj. Therefore, in Eq. (j20J, 
approximating P by a s P SL ^ should lead to a better ap- 
proximation for the evolution than approximating P by 
zero. Not surprisingly, with this approximation, Eq. I|29(l 
can be regarded as a generalized version of the MLLA 
equation to include quarks, in the sense that the g com- 
ponent of this latter equation for D = (P/£,D S ) when 
Eq. I|47|) is invoked is precisely the MLLA equation of 
Ref. |3 We therefore conclude that although Eq. lj29*|) 
is derived from the DLA, it is more complete than the 
MLLA equation since in Eq. 129|) it is neither necessary 
to restrict P in this way nor to use the approximation in 
Eq. g7J). 

We may now approximately but explicitly calculate 
the evolution in the SGL+FO(+FO<5) scheme by approx- 
imating P SGL in Eq. I|14l) by its leading term P DL given 
by Eq. O, being the first, 0(a s /uj), term in Eq. (|H)j) . 
Thus we take 



P = P 



DL 



(51) 



Recall that P FO is equal to P in Eq. (JUJ) when all SGLs 
are excluded. After inverse Mellin transforming Eq. H51fl . 
we can solve Eq. Q directly using standard numerical 
techniques for x space evolution, as suggested in Sec. 
Hfl We shall call this the DL+FO(+FO<5) scheme. We 
note that, although the analytic solution to Eq. @ for 
P = P DL as defined in Eq. g2> is given by Eq. I|44|l . 
while the analytic solution to Eq. © for P = P FO can be 
found using the well known method in the FO approach, 
an analytic solution to Eq. © for P given in Eq. (|5 1|) 
does not seem feasible. 

For the DL+FO(+FO<5) scheme, we require 
P Dh (x,a s ). The inverse Mellin transform of Eq. 
O gives 

P DL (x,a s ) = ( 4 v ^7lni') , (52) 

x m - V x J 

where Ji(y) is the Bcsscl function of the first kind, given 
by 



t ( \ V V 4 / 

, — n \ ' 



(53) 



This gives 



(~l) r 
r!(r + l)! 



(4C> S In 2 x) r 



(54) 

The series in Eq. (|54l) may also be obtained by expanding 
Eq. 142(1 to infinite order in a s , using Eq. I|13|) to perform 
the inverse Mellin transform on each term, and finally 
using the identity 

( ~ 4) J2r)\ " H (55) 

for r > 0. Equation (|54|l in fact converges rapidly, how- 
ever the truncated series can differ substantially from the 
full series whenever x is small enough. For more reliabil- 
ity, the Bessel function should be calculated by numerical 
evaluation of the result 



My) = - I M cos( y sine -6). 
k Jo 



(56) 



Note that for Q > Aqcd, the evolution of D(x,Q 2 ) in 
the DL+FO(+FO<5) scheme contains no Landau pole for 
all x, not just for x > Aqcd /Q as is the case if we evolve 
with Eq. l(25)l. 

Using the property at large y that 



Jl<!/) = V^ co8 ( ! '"T I 



(57) 
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we find the true x — » divergence 



P UL (x,a s ) 



(C A a s )*A 



+ o 




cos 



4^ C^fts In - - V 
a; 4 



(58) 



one from the FO 
is proportional to 



of P to be weaker than the 
approach, which from Eq. (|54|l 

(l/x)In 2(n-1) (l/aj) at 0«). 

For simplicity, in the next section we shall take P FO 
(and P*2q) to 0(a s ) only in the DL+FO(+FO<5) scheme. 
In other words, we will approximate P by 



P = P DL + a s P F °^(+a s P 



SL(0h 



(59) 



We shall call this the DL+LO(+L05) scheme. 

It is well known from the MLLA that while DLs give 
the shape of the peak that occurs at small x, the gluon 
component of the type p = SL is required to get the 
correct peak position. We therefore anticipate that, rel- 
ative to both the DL+LO scheme and the MLLA, the 
DL+LO+L05 scheme will give a better description of 
the data since it contains both gluon and quark compo- 
nents of the type p = SL. 

In Fig.^ we see that P gg (x, a s ) in the DL+LO scheme, 
which is equal to the DL+LO+LO(5 scheme when i^l, 
interpolates well between its 0(a s ) approximation in the 
FO approach at large x and P® g h (x,a s ) at small x (the 

small difference here comes from P FO (°)(:x;) at small x). 
DL resummation clearly makes a large difference to P at 
small x. 



IV. COMPARISONS WITH DATA 

In this section, we elaborate on our numerical study 
of the DL+LO+LO(5 scheme in Ref. |Io|. We perform 
a numerical comparison of the FO and DL+LO+LO<5 
schemes by otherwise imposing the same assumptions 
and choice of parameterization, starting scale etc., and 
then fitting in each scheme to precisely the same experi- 
mental measurements of the normalised differential cross 
section for light charged hadron production in the pro- 
cess e + e~ — > (7, Z) — > h + X, where h is the observed 
hadron and X is anything else. This data, spread over a 
wide range in center-of-mass energy y/s, is com pos ed of 
the sets from TASSO at j£ = 14, 35, 44 GeV and 
22 GeV MARK II and TPC at 29 GeV, 
TOPAZ at 58 GeV [13, ALEPH [H|, DELPHI 0, L3 
[il, OPALJ21I and SLC at 91 GeV, ALEPH E| 
and OPAL || at 133 GeV, DELPHI at 161 GeV El 
and OPAL at 172, 183, 189 GeV ;26| and 202 GeV [U 
These data span a wide range in x p = 2p/ -^/s, where p is 
the momentum of the observed hadron, which constrain 



- DL+LO(+L05) 

DL 
- LO 



5 
4 
3 
2 

1 W 
- 
^-1 - 
-2 - 
-3 - 
-4 - 
-5 - 

-6 I ' ^ 

01 234567 

In (1/x) 

FIG. 1: (i) P gg (x,a s ) calculated in the DL+LO(+L05) 
scheme, (ii) P gg (x,a 3 ) calculated to 0(a a ) in the FO ap- 
proach (labelled "LO"), and (iii) P 99 L (:r,a s ) (labelled "DL"). 
a s = 0.118/(2tt). 



the FFs in the region of x for which x p < x < 1. In Ref. 
[l2|. large £ data could be described within the range 



£<ln 



2M' 



(60) 



where £ = ln(l/a;j,) and M is a mass scale of 0(1) GeV. In 
Ref. Q, this inability to describe large £ (small x p ) data 
using the naive approach to the MLLA was attributed 
to the formally beyond-MLLA evolutions of the higher 
moments (skewness, kurtosis etc.) becoming too large at 
large a/s, and when these were fixed to be zero it was 
found that the MLLA could give an excellent description 
of data up to the highest values of £ currently measured. 
In the DL+LO+LOJ scheme here, the evolution of the 
higher moments should be somewhat suppressed by the 
FO contribution [9|. However, we will nevertheless im- 
pose Eq. I|60l) on the data to be fitted to, since here our 
aim is to extend the good FO DGLAP description of 
small £ (large x p ) data to larger £ via DL resummation. 

At LO in the coefficient functions, these data are de- 
scribed in terms of the evolved FFs by 



1 da 
a(s) dx p 



1 



n f (Q(s)) 



Q q (s)D+(x p ,Q 2 ), (61) 



where Q q is the electroweak charge of a quark with 
flavour q and (Q) is the average charge over all flavours. 
For Eq. (|61|l to be a valid approximation, it is necessary 
to choose Q = 0(y/s). Since we only use data for which 
V~s > nif,, where mt, « 5 GeV is the mass of the bot- 
tom quark, and since later we will also set Qq > nib, 
we will take rif = 5 in all our calculations. While the 
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precise choice for n / does not matter in the DLA, calcu- 
lations in the FO approach depend strongly on it. Since 
we sum over hadron charges, we set Dq = D q . Since we 
do not use data with quark tagging, the c quark cannot 
be distinguished from the u quark since both quarks cou- 
ple to the Z boson in the same way, i.e. have the same 
electroweak charge. Likewise, the d, s and b quarks are 
similar to one another in this respect. Therefore, to avoid 
redundant degrees of freedom, we fit only the FFs 



fuc(x, Qo) = 2 i u ( x > Qo) + c ( x > Ql)) > 
fdsb(x, Ql) = \ {d(x, Ql) + s(x, Ql) + b(x, Ql)) 



(62) 



and the gluon g(x,Q 2 ) ). For each of these three FFs, we 
choose the parameterization 



f(x,Ql) = N exp[-c In 2 x]x a (1 - x) 



13 



(63) 



since at intermediate and large x the FF is constrained 
to behave like 



f(x,Q 2 )^Nx a (l-xf, 



(64) 



which is the standard parameterization used in global fits 
at large x, while at small x (where (1 — x)@ ss 1) the FF 
is constrained to behave like 



lim f(x,Q ) = iVexp 



1 2 1 , 1 

-cln a in — 

x x 



(65) 



which for c > is a Gaussian in ln(l/x) of width 
l/\/2c, center at —a/ (2c) and normalization given by 
iV \Ztt I cexp [a 2 /(4c)] . For sufficiently large Qo, the DLA 
predicts such behaviour with a < 0. 

In addition, we use Eq. to remove four free pa- 
rameters by imposing the constraints 



Cue " 



Cdsb = 
z Oidsb 



Cg, 
= OL a 



(66) 



This implies that all FFs have the same width and cen- 
ter, however the normalisations may not turn out to be 
consistent with Eq. I|48|l . However, Eq. I(48|l is only an 
approximation at small x, while the N are also relevant 
in the large x region, where also the (1 — x) 13 factors are 
necessary. At any rate, it will be interesting to see how 
well the relation 



Cw 

N uc « N dsb ^—N g 



(67) 



as implied by Eq. 1)48(1 is obeyed after a fit is performed. 
In addition to the 8 free parameters for the FFs, we also 
fit Aqcd- We choose Q 2 — s, although it is only im- 
portant that the latter two quantities are kept propor- 
tional, since the constant of proportionality has no effect 
on the final FF parameters and the description of the 
data (or, equivalently, the quality of the fit). However, 



the final fitted Aqcd varies in proportion to this con- 
stant, so there will be an overall theoretical error on our 
fitted values for Aqcd of a factor of O(l). Since all data 
will be at t/s > 14 GeV, we choose Qo = 14 GeV. As 
discussed in Sec.[H] the evolution is performed by numeri- 
cally integrating Eq. ■ For this we use a grid consisting 
of 250 points equally spaced in InQ 2 over the range 14 
GeV< Q < 202 GeV, and 750 points equally spaced in 
ln(l/a;) over the range < \n(l/x) < 11.6. 



A. Fixed Order Evolution 

We first perform a fit to all data sets listed above using 
DGLAP evolution in the FO approach to LO, without DL 
resummation. This approach is the same as that used in 
fits in the literature. We fit to those data for which Eq. 
(|rJ0j) is obeyed with M = 0.5 GeV. This gives a total of 
425 data points out of the available 492. We obtain a x 2 
per degree of freedom Xdf = 3-0 (or 2.1 after subtraction 
of the contribution to \ 2 from the TOPAZ data, which is 
the only data set from which an individual Xbf greater 
than 6 is obtained), and the results are shown in Fig. 
121 and Table [I] The result for Aqcd is quite consistent 
with that of other analyses, at least within the theoretical 
error. It is clear that FO DGLAP evolution fails in the 
description of the peak region and shows a different trend 
outside the fit range. The exp[— cln 2 x] factor does at 
least allow for the fit range to be extended to x values 
below that of x — 0.1, the lower limit of most global fits, 
to around x = 0.05 (£ = 3) for data at the larger 
values. Note that the negative value of j3 for the gluon is 
unphysical, because the gluon FF is weakly constrained 
in our fit since it couples to the data only through the 
evolution (see Eq. 1)61(1 '). 



TABLE I: Parameter values for the FFs at Q = 14 GeV 
parameterized as in Eq. (1631 from a fit to all data listed in 
the text using DGLAP evolution in the FO approach to LO. 
Aqcd = 388 MeV. 



~~~~~--~--__I 5 ar ameter 


N 


P 






FF 


o 


c 


9 


0.22 


-0.43 


— 2.38 


0.25 


u + c 


0.49 


2.30 


[-2.38] 


[0.25] 


d + s + b 


0.37 


1.49 


[-2.38] 


[0.25] 



B. Incorporation of Soft Gluon Resummation 

We now perform the same fit, i.e. to the same data 
with the same parameterization, but now evolving in the 
DL+LO+LO<5 scheme. The results are shown in Table 
|H|and Fig. □ We obtain Xdf = 2A ( or lA without 
the TOPAZ data — for each remaining data set, the in- 
dividual Xbf i s ^ ess than 3), a significant improvement 
to the fit above with FO DGLAP evolution. The data 
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FIG. 2: Fit to data as described in Table|I] Some of the data 
sets used for the fit are shown, together with their theoretical 
predictions from the results of the fit. Data to the right of 
the horizontal dotted lines have not been used in the fit. Each 
curve is shifted up by 0.8 for clarity. 



TABLE II: Parameter values for the FFs at Q = 14 GeV 
parameterized as in Eq. I)6M|I from a fit to all data listed in the 
text using DGLAP evolution in the DL+LO+L05 scheme. 
Aqcd = 801 MeV. 



~~~~~--~--__I 5 arameter 


N 


P 






FF 


cx 


c 


9 


1.60 


5.01 


— 2.63 


0.35 


u + c 


0.39 


1.46 


[-2.63] 


[0.35] 


d + s + b 


0.34 


1.49 


[-2.63] 


[0.35] 



around the peak is now much better described. The en- 
ergy dependence is well reproduced up to the largest yfs 
value, Vs = 202 GeV. At ^fs = 14 GeV, the low x de- 
scription of the data is extended from x = 0.1 in the un- 
resummed case down to 0.06 in the resummed case, and 
from x = 0.05 to 0.005 at y/s = 202 GeV. This should 
also be compared to the fit to the same data in Ref. 
[l2| where DL resummation was used within the MLLA 
but with neither FO terms nor quark freedom (i.e. Eq. 
(I48|l was imposed over the whole x range) . That fit gave 
X DF = 4.0. We conclude that, relative to the MLLA, the 
FO contributions in the evolution, together with freedom 
from the constraint of Eq. (|48|) . makes a significant im- 
provement to the description of the data for £ from zero 
to just beyond the peak. 

The value Aqcd ~ 800 MeV is somewhat larger than 
the value 480 MeV which we obtain from a DGLAP fit in 
the large x range (x > 0.1). We note that had we made 
the usual DLA (MLLA) choice Q = \fs /2 instead of our 
choice Q — \fs which is usually employed in analyses 
using the DGLAP equation, we would have obtained half 



FIG. 3: Fit to data as described in Table El 



this value for Aqcd . N g is too large by a factor of about 
2 relative to its prediction in Eq. i|67|). As noted before, 
the initial gluon FF is weakly constrained in our fits since 
it only enters the cross section via the mixing with the 
quark FFs in the evolution. 

While the DL resummation greatly improves the de- 
scription around the peak, the description still deteri- 
orates as £ increases, since a fit with data at larger £ 
resulted in an increase in Aqcd and N g , as well as Xdf- 
The large £ data may be better described with the in- 
clusion of the (unknown) resummed SL contribution. In- 
deed, fitting without the LO p = term (the LO SL), i.e. 
fitting in the DL+LO scheme, generally gives a larger 
Xq F compared to the same fits in the DL+LO+L05 
scheme, since the data around the peak region could not 
be described. 

We repeated the above two fits without the constraints 
given in Eq. (|66[) and found no substantial improvement 
to the fits. 

The results of these two fits show that, for the DGLAP 
evolution to describe the overall features of the data from 
small to large x p , resummation of SGLs is necessary. In 
addition, our DL resummation scheme proves to be an 
adequate implementation of this resummation. However, 
the description at £ values beyond the peak remain in- 
accessible. It is clear from the smoothness of the shape 
around and beyond the peak in Fig. [3] that this cannot 
be attributed to any instability in the evolution of the 
higher moments. This instability occurred in the method 
of Ref. [12, which was pointed out and remedied in Ref. 
|9J. Figure in fact shows that the evolution at large £ 
follows the data well. Indeed, a fit in the DL+LO+LO<5 
scheme to data in the same region as that used in the fits 
of Ref. 0, 



£ > 0.75 + 0.33 ln(Vs), 



(68) 



10 



gives Xdf = 7, which is admittedly unacceptable large, 
however an FO fit to the same data gives Xdf = 35- Thus 
there is a problem with simultaneously describing both 
the range from £ = to the peak region to the range be- 
yond the peak within the DL+LO+LCW scheme, imply- 
ing a large theoretical error beyond the peak. This may 
be due to the neglect of the complete SL contribution, 
which is unfortunately unknown. Another likely reason 
is our neglect of the effect of the produced hadron's mass. 
This effect is important at small x p . We will therefore 
study this effect in the next subsection. 



C. Incorporation of Hadron Mass Effects 

We will now incorporate hadron mass effects into our 
calculations, using a specific choice of sca ling variable. 
An alternative approach is given in Ref. [28( (see also 
Ref. |29j). For this purpose it is helpful to work with 
light cone coordinates, in which any 4- vector V is written 
in the form V = (V+, V~ , V T ) with V ± = -^{V° ± V 3 ) 

and Vt = (V 1 , V 2 ). In the center-of-mass (COM) frame, 
the momentum of the electroweak boson takes the form 







(69) 



In the absence of hadron mass, x p (whose definition 
x.p = Ipj^fs applies only in the COM frame) is iden- 
tical to the light cone scaling variable r\ = pt How- 
ever, the definition x p = 2p/^/s applies only in the COM 
frame, so 77 is a more convenient scaling variable for 
studying hadron mass effects since it is invariant with 
respect to boosts along the direction of the hadron's spa- 
tial momentum. Taking this direction to be the 3-axis, 
and introducing a mass for the hadron, the momen- 
tum of the hadron in the COM frame reads 



Ph 



(70) 



Therefore the relation between the two scaling variables 
in the presence of hadron mass is 



V 



(71) 



Note that these two variables are approximately equal 
when m,h <C x p \fs, i.e. hadron mass effects cannot be 
neglected when x p (or rf) are too small. 

In the leading twist component of the cross section af- 
ter factorization, the hadron h is produced by fragmen- 
tation from a real, massless parton of momentum 



k 



ElL 

y 



0,0 



(72) 



The + component of everything other than this parton 
and of everything produced by the parton other than 



the observed hadron h must be positive, implying y > 
77 and y < 1 respectively. As a generalization of the 
massless case, we assume the cross section we have been 
calculating is (da/dr))(ri,s), i.e. 



da 
drj 



1 dy da 
V dy 



(y,s,Q 2 )D(^,Q 



(73) 



which is related to the measured observable 
{da / dxp){x p , s) via 



da 



1 



da 



1 



sri 2 (x p ) 



drj 



( V (x p ),s). (74) 



Note that the effect of hadron mass is to reduce the 
size of the cross section at small x p (or 77), which Fig. 
suggests is what is needed to improve the fit. The 
problem with the above method is that it only applies 
to the case where only one species of hadron is pro- 
duced, whereas the data we are studying are for light 
charged hadrons, i.e. charged pions, charged kaons and 
the (anti)protons, whose masses (140, 494 and 938 MeV, 
respectively) are substantially different. However, since 
most of the produced particles will be pions and kaons, 
the effective range of hadron masses may be sufficiently 
small to achieve reliable results when the particle masses 
are taken to be equal, which we will do. We now per- 
form the last two fits again but with mh, included in 
the list of free parameters. For FO evolution, we obtain 
Xq F = 2.06, which is a substantial improvement over the 
same fit above for which no treatment of hadron mass 
effects is applied. The results are shown in Fig. Q] and 
Table ITTT1 The result for m/j is of the expected order of 
magnitude, however Aqcd is unreasonably large. The 
suppression of the cross section beyond the peak from 
hadron mass effects is evident, and allows for the cross 
section to follow the data much more closely. 



TABLE III: As in Table [I] but incorporating mass effects in 
the fit. Aqcd = 1308 MeV and m h = 408 MeV. 
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FF 
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0.11 


-0.82 


— 2.01 


0.18 


u + c 


0.70 


2.12 


[-2.01] 


[0.18] 


d+ s + b 


0.82 


2.35 


[-2.01] 


[0.18] 



For the DL+LO+LO<5 fit, we obtain the results in Fig. 
and Table IIVI The parameters are not substantially 
different to those in Table ITT1 The result for mt is again 
reasonable. We find Xdf = 2.03, i.e. the quality of the fit 
is the same as for the previous FO fit, showing that the 
FO case benefits much more from the inclusion of mass 
effects than the DL resummed case. However, treatment 
of mass effects renders the value of Aqcd obtained in the 
fit with DL resummation more reasonable. 

We conclude therefore that to improve the large £ de- 
scription and to achieve a reasonable value for Aqcd, 
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resummed DL contribution given by Eq. (|52|l must be 
added to the NLO splitting functions. These contain 
SGLs belonging to the classes m = 1, 4, which must be 
subtracted. Note that the NLO m = 1 term is accounted 
for by the resummed DL contribution. The m = 4 term 
is a type p = term, and hence does not need to be 
subtracted. 



D. Comparison With Gluon Jet Data 

We now compare our results with the OPAL gluon jet 
measurements at E- )ct = 14.24, 17.72 [30j and 40.1 GeV 
[3l| . These data are shown in Fig. [U together with our 
gluon FF from the fit of Table [H] incorporating DL re- 
summation. The initial gluon FF from that fit is found 
to be about twice as large as predicted by Eq. I|67|l . and 
this is reflected in the figure. In Fig. [7] we show the same 



FIG. 4: Comparison with data from the fit described in Table 

ITTT1 



TABLE IV: As in Table ITT1 but incorporating mass effects in 
the fit. Aqcd = 399 MeV and m h = 252 MeV. 
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FF 




c 
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1.59 


7.80 


— 2.65 


0.33 


u + c 


0.62 


1.43 


[-2.65] 


[0.33] 


d + s + b 


0.74 


1.60 


[-2.65] 


[0.33] 




FIG. 5: Comparison with data from the fit described in Table 

hyi 



both DL resummation and treatment of mass effects are 
required. 

To extend NLO calculations to small x p , the complete 




FIG. 6: Comparison of the gluon FF in the DL+LO+LOS 
scheme to gluon jet measurements from OPAL. The postion 
of the cut of Eq. 16011 is shown to indicate where the gluon is 
constrained. Q = 2_Ej ct . 

plot again, but this time from a fit in which these data 
are included (but which is otherwise identical to the fit of 
Table InJ. The data is identified with our evolved gluon 
FF at Q = 2_Ej Ct . In this case Xdf = 2-3, a slightly larger 
value than that from the fit of Table ITT1 in particular be- 
cause the gluon jet data at Ej Ct = 14.24 GeV cannot not 
be well fitted. These data give an individual Xdf °f 3-7, 
although the data at E iet = 17.72 and 40.1 GeV give 0.9 
and 1.2 respectively. The data around and beyond the 
peak are poorly described. The parameters and compar- 
ison with the remaining data are shown in Table and 
Fig. [S] The value for N g is lower than that in Table 
ITT1 and in better agreement with Eq. 1)67(1. This smaller 
gluon FF is presumably the cause of the undershoot of 
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the calculation from the data at high y/s seen in Fig. [S] 
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FIG. 8: Fit to data as described in Table [V] 



FIG. 7: As in Fig. but from a fit in which the gluon jet 
data below the cut are included in the fit. 



TABLE V: Parameter values for the FFs at Qo = 14 GeV 
parameterized as in Eq. 163H from a fit to all data, includ- 
ing the OPAL gluon jet data using DGLAP evolution in the 
DL+LO+LOS scheme. A QC d = 954 MeV. 
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1.27 


4.21 


— 2.41 


0.29 


u + c 


0.40 


1.55 


[-2.41] 


[0.29] 


d + s + b 


0.50 


1.59 


[-2.41] 


[0.29] 



Repeating the fit with hadron mass effects accounted 
for (including in the gluon jet data), we obtain Aqcd = 
490 MeV, m h = 302 MeV, and Xdf = 2 - L The fitted 
parameters, shown in Table IvTl are not significantly dif- 
ferent to those of TableEl The comparison with the data 
fitted to are shown in Figs. 151 and 1 101 The individual Xbf 
values for the gluon jet data at E- ]ct — 14.24, 17.72 and 
40.1 GeV are now 1.2, 0.4 and 0.7 respectively. We con- 
clude that the description of the gluon jet data are also 
improved by including hadron mass effects, particularly 
at low -Ejet- 

In Fig. 03 we observe that the experimental data have 
a tendency to systematically undershoot predctions at 
large £. This could be partially caused by the different 
procedures for including soft particles in the definitions of 
quark and gluon jets used in the experimental analysis. 
We observe furthermore that these deviations decrease 
with increasing energy. 




FIG. 9: As in Fig. but with a fitted hadron mass. 



TABLE VI: Parameter values for the FFs at Q = 14 GeV 
parameterized as in Eq. I|63fl from a fit to all data, includ- 
ing the OPAL gluon jet data using DGLAP evolution in the 
DL+LO+L05 scheme and with mass effects incorporated. 
Aqcd = 490 MeV. 



~~~~-------_I 3 arameter 


N 


P 






FF 


o 


c 


9 


1.30 


5.09 


— 2.30 


0.24 


u + c 


0.46 


1.70 


[-2.30] 


[0.24] 


d+s + b 


0.53 


1.75 


[-2.30] 


[0.24] 



13 



20 



15 - 



♦ OPAL 202 
i OPAL 1 72 

a DELPHI 161 

* OPAL 133 
. OPAL 91 

□ TASSO 44 

* TASSO 35 

• TPC 29 
TASSO 14 



■3 10 - 




FIG. 10: Fit to data as described in Table IVTI 



V. CONCLUSIONS 

We have defined a general scheme for resumming SGLs 
in DGLAP evolution to any order in the FO contribu- 
tion and to any class in the resummed SGL contribu- 
tion. We have implemented it numerically at LO in the 
FO contribution with DL resummation, using explicit re- 
sults on the LO splitting functions from the literature. 
This DL+LO+L06 scheme is obtained by taking the LO 
results for the splitting functions, subtracting the un- 
resummed DLs and adding the complete resummed DL 
contribution, given by Eq. (|52|) , which we obtained using 
the DLA. This scheme contains all contributions from the 
MLLA, as well as all other FO contributions at LO. We 
have shown by fitting to data with this scheme that for 
DGLAP evolution to also describe large £ data around 
the peak region, resummation of SGLs is necessary. In 
addition, we showed that both SGL resummation and 
treatment of hadron mass effects leads to a good fit and 
a reasonable result for Aqcd- 

In Ref . [l^l , it was shown that the MLLA alone allows 
for a good description of data around the peak region. 
In Ref. Q, it was shown that provided certain spurious 
higher order terms were removed, data above the peak 
region could also be well described. However, the low 
£ region was not well reproduced, while in Ref. |12( the 
good description here was due to a well known coinci- 
dence. Simultaneously fitting to both small and large 
£ data using either of the two approaches leads to ex- 
tremely high Xt>f values. Using the approach in this 
paper gives a much better fit to all the data, even if the 
fit is still not in the acceptable range. It allows for the 
data in a larger £ range to be described than the FO ap- 
proach and the MLLA approach of Ref. 01 do. Further 
improvement in the large £ region can be expected from 



the inclusion of higher order SGLs. 

Our scheme allows a determination of quark and gluon 
FFs over a wider range of data than previously achieved, 
and should be incorporated in NLO global fits of FFs 
such as that in Ref. |2j by using the method of the last 
paragraph of Sec. IIVI since the current range of 0.1 < 
x p < 1 is very limited. 



APPENDIX A: TIMELIKE SPLITTING 
FUNCTIONS IN PERTURB ATIVE QCD 

The LO terms of the timelike splitting functions in x 
space for the component D = (-Ds, D g ) are 



P^(x) =C F l-l-x + 2 



1 



1-x 



+ ~£(l-x)J, 



(0) 



2C, 



1 + (1 - xf 
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(0) 



p(0) 
39 



T R n f (x 2 + (l-x) 2 ), 



2C A - - 2 + x - x 1 
\ x 



1-x 



+ (yC A - |?W) 6(1 - x), 



(Al) 



where Tr = 1/2, rif is the number of flavours and, for 
the color gauge group SU(3), Ca — 3 and Cf = 4/3. 
The function [1/(1 — x)] + is defined by 



dy 



1-2/ 



/(- 



11 
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1-1/ 



/ ( - ) - yf{x) 



f{x) ln(l - x) 
(A2) 

for any function f(x). Transforming Eq. (|A1|) to Mellin 
space gives 



1 



P^) = T R n f 
P^( U ) = 2C A 



2 (w + l)(w + 2) 

J 1 + 3w + 4 
lu(uj + 1)(lu + 2)' 
u? + 3oj + 4 
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-Si(u> + 1) 
where, for integer n, 



n , 



(A3) 
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fe=i 
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The DLs and SLs at LO are obtained by expanding 
Eq. (IA3|) about uj = 0, for which the result 

S 1 (w + l) = l + 0(u}) (A5) 

is required. 

The NLO splitting functions P^(x) are presented in 
Ref. j^, while their Mellin transforms pW(w) are pre- 
sented in Ref. 0. 

We do not explicitly present p( 0,1 \x) for the compo- 
nents D = D~ , Dns, since it is enough for our purposes 



to know that they do not contain SGLs. 
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